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COMMUNICATION CHANNEL CAPACITY ESTIMATION 

Related Application Data 
[0001] This application claims the benefit of and priority under 35 U.S.C. § 1 1 9(e) to U.S. 
Provisional Application No.: 60/543,966, entitled "Equivalent Estimation Method for 
Evaluating Subscriber Lines Based On Time Domain Reflectometry," filed February 1 1, 
2004, which is incorporated herein by reference in its entirety. 

BACKGROUND 

Field of the Invention 

[0002] This invention generally relates to communication systems. In particular, an 
exemplary embodiment of this invention relates to estimating communication chaimel 
capacity. More particularly, an exemplary aspect of the present invention relates to 
estimating digital cormnunications channel capacity. 

Description of Related Art 

[0003] Digital Subscriber Line (DSL) technology makes it possible to transport high bit 
rate digital information via a commxmications channel, such as a subscriber line. The channel 
capacity, which is defined as the obtainable data rate for a given line, is based on the physical 
structure and topology of the line, such as the length of the line, gauge, existence of bridged 
taps, bridged tap locations and lengths, etc.. Thus, if the topology of the line is known, the 
data rate can be predicted prior to providing DSL service to a customer, 

J 

[0004] Time domain reflectometry (TDK) is a very useful tool for characterizing a 
subscriber line. TDR operates by sending an electrical pulse down the line and measuring the 
returned signal, referred to as a TDR echo. The measured TDR echo contains information 
about the physical structure and topology of the line. 

[0005] The most common method for evaluating the channel capacity from a TDR echo 
is by explicitly estimating the physical topology of the line based on the transmission line 
theory, and then searching a database to find a data rate corresponding to the specific 
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topology. One issue with this approach is the estimation complexity increases dramatically 
when the topology of the line is complex. For example, consider a line with N consecutive 
sections of dilBferent gauges but without bridged taps. Since there are two variables that need 
to be estimated for each section, the gauge and the length, the searching space is 2N 
dimensional. Thus, the computation complexity increases exponentially with the nimiber of 
distinct sections. 

SUMMARY 

[0006] The ultimate goal of estimating the topology of a subscriber line from the TDR 
echo is to evaluate the channel capacity of the line. Hence, if one can estimate a line with a 
much simpler topology but the same channel capacity as the actual line from the TDR echo, 
the goal can be achieved with much less computation complexity. 

[0007] An exemplary embodiment of the present invention follows this theory, and 
provides a method of evaluating the channel capacity from the TDR echo without explicitly 
estimating the topology of the line. 

[0008] For simplicity, herein a subscriber line without bridged tap is referred to as a 
straight line. Here "straight" means no bridged tap presented. Also, a subscriber line is 
referred to herein as a line, a loop or a wire interchangeably. 

[0009] An exemplary embodiment of the present invention is developed based on 
transmission line theory. A typical transmission line system can be schematically represented 
as illustrated in Fig. 1, where Vs is the output voltage of the source, Zs is the output impedance 
of the source, Zl is the input impedance of the load, and L is the length of the line connecting 
the source and the load. Vs, Zs and Zi are functions of frequency. According to the 
transmission line theory, when the electrical wave generated by the source travels down to 
the line, it is attenuated by the line, and reflected whenever there is an impedance 
discontinuity. The wave travels forward and backward inside the line to infinity. 

[0010] Assume the characteristic impedance and the propagation constant of a given line 
to be Zo and where Zq and j^aie functions of frequency. With a straight loop with a single 
gauge, the only impedance changes are (1) at the cozmection between the source and the line, 
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and (2) at the connection between the line and the load. Assume the reflection coefficient, 
which is defined as the ratio between the backward voltage wave (toward the source) and the 
forward voltage wave (away firom the source), to be at the source output, and to be pL at the 
load, where 



Pl -~ 



[0011] According to the wave propagation theory, the voltage at the output of the source, 
denoted as F©, can be represented by a combination of an incident wave and infinite number 
of multi-reflections caused by impedance discontinuity: 



Incident wave: 



1^ reflection: 



2"^ reflection: 



3"* reflection: 



^0+2-, 



Forward wave V^^^ = Vj°^ • e"^*^ ' Pl' Ps y 
Backward wave V^^ = v}'^^ • e'^^ ■ pj^ , 

Forward wave V^^^ = Fg^") • e^"^ • p\ - p] , 
Backward wave V^^>- = V}°^ • e"^'*' • p\ ' p^ , 

Forward wave V^^'^ = • e"*'*' • p\ • p\ , 
Backward wave V^""^ = fJ") • e^^ • • p] , 



n* reflection: Forward wave V^"^ = V^"^ • e'^'^ • p"^ • p^ , 

Backward wave Vj"^ = V}"^ ■ e'^"* • p"^ ■ p^-^ . 
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[0012] 



Therefore, Vo can be e^spressed as: 



(1) 



1 + e 



,-2yi 



Pl 



[0013] Likewise, the voltage at the input of the load, denoted as Vl, can also be 
represented by a combination of multi-reflections: 



1* reflection: Forward wave J^'^ = v}"^ • e"'' , 

Backward wave = • e"'' • p^ , 

Forward wave vj;^^ = v}''^ • e-^^ Pl'Ps > 
Backward wave V^^^ = fJ") • e"^'* ' p\- Ps ^ 

Forward wave V^^ = Fo^") • e"*'^ • /? J • p] , 
Backward wave Fp>- = v}''^ • e"^'^ • P^ • p] , 



2nd 

reflection: 



reflection: 



n*^ Forward wave V^f^ = Fo^»> • e'^^'-'^^ • • , 

reflection: ' 
Backward wave V^^ = pKo) . g-2(«-i);i .^/. 



Therefore, Vt can be expressed as: 
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[0014] Both a TDR measurement system and a DSL application can be represented by 
flie equivalent circuitry shown in Fig. L In the TDR case, Vs is Ihe pulse sent down to the 
line, Vo is the measured TDR echo. The iucident wave in Vq is referred to as the near-end 
echo, the sum of the multi-reflections is referred to as the far-end echo. In a TDR 
measurement, the source impedance is usually the same as the characteristic impedance of the 
line, i.e. Zs = Zq, and the end of the line is usually open, i.e. Zi = a?, thus Oy pl — 1^ and 
the measured TDR echo can be written as: 

^o=^."7^-{l + ^-''') (3) 
which is a combination of the near-end echo and the backward wave in the 1^* reflection. 

[0015] In the DSL application case, Vs and Zs represent the eqxiivalent circuitry of the 
modem at the central office (CO), and Zl represents the equivalent circuitry of the modem at 
the customer premise (CPE). The obtainable data rate relates to the transfer function of the 
subscriber line, which is defined as: 

K 

[0016] A modotn is usually designed to have an impedance matching to the line, i.e. Zs = 
Z(h Zl = Zo, thus Ps = 0, pl = 0, and Vo, Vl, and H can be written as: 

V„ =V — ^2 

V^=K — e-^, 

H^e-'^. (4) 

[0017] Because the imaginary part of is a linear function of frequency, H has linear 
phase, the data rate is mainly determined by the modular of the transfer function. 

[0018] Assume tiie transfer functions of two single-gauge straight lines to be Hj and H2. 
According to Eq. (4), 
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where and Li are the propagation constant and length of Line 1, respectively, and y2 and L2 
are the propagation constant and length of Line 2, respectively. These two lines will have the 
same data rate if: 

realiy^)' = real{y^)Lj^ , (5) 
where real(') is an operation to obtain a variable's real part 

[0019] Assume the TDR echoes of these two lines to be Vqj and ¥02- According to Eq. 
(3), 

where Zoi and Z02 are the characteristic impedance of Line 1 and Line 2. Ignoring the 
difference in characteristic impedance between these two lines, i.e. Zoj ^Zo2 - Zo> these two 
lines have sunilar near-end echoes, then according to Eq. (5), when these two lines have the 
same data rate, the following equation holds true: 



'Z.-hZ. 



(6) 



[0020] The left-hand side of Eq. (6) is the amplitude of the far-end TDR echo of Line 1 ; 
the right hand side is that of Line 2. Eq. (6) mdicates that the far-end TDR echoes from two 
single-gauge straight lines, which have the same data rate, and have the same amplitude. 
Although this deduction is derived for certain values of Zs and Z/,, the conclusion is 
applicable for any Zs and Zl. 



[0021] An exemplary embodiment of the present invention is developed based on this 
deduction. 



6 



wo 2005/078951 PCTAJS2005/004015 
BRIEF DESCRIPTTON OF TRE DRAWINGS 



[0022] Fig. 1 illustrates a schematic representation of a transmission line system with a 
single-section strai^t loop; 

[0023] Fig. 2 illustrates the exemplary spectrum of the TDR pulse used according to an 
exemplary embodiment of this invention; 

[0024] Fig. 3 illustrates exemplary spectra of the far-end echoes from 26AWG straight 
loops with various lengths according to this invention; 

[0025] Fig.4a illustrates a data rate comparison between equivalent-data-rate loops for the 
US data rate according to this invention; 

[0026] Fig.4b illustrates a data rate comparison between equivalent-data-rate loops for 
the DS data rate according to this invention; 

[0027] Figure 5a illustrates the far-end echo comparison between 24AWG and 26AWG 
equivalent-TDR-echo loops where the 26AWG equivalent loop length = 5kft. according to 
this invention; 

[0028] Figure 5b illustrates the far-end echo comparison between 24AWG and 26AWG 
equivalent-TDR-echo loops where the 26AWG equivalent loop length = lOkfl. according to 
this invention; 

[0029] Figure 6 illustrates schematically a representation of a transmission line system 
with a 2-section straight loop; 

[0030] Figures 7a and 7b are a far-end echo comparison between mixed-gauge loops and 
26AWG equivalent loops according to this invention; 
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[0031] Figures 8a and 8b illustrate a data rate comparison between 2-section mixed- 
gauge loops and their 26AWG equivalent-data-rate loops for the US data rate and DS data 
rate, respectively, according to this invention; 

[0032] Figures 9a and 9b illustrate the data rate comparison between 26AWG+24AWG 
loops and 24AWGrf 26AWG loops for the US data rate and the DS date rate, respectively, 
according to this invention; 

[0033] Figurel 0 illustrates a straight loop with mixed gauges according to this invention; 

[0034] Fig. 1 1 illustrates an exemplary channel capacity estimator 100 according to this 
invention; 

[0035] Figures 12a and 12b illustrate the loop length estimation results, with Fig. 12a 
showing a comparison between the actual and the estimated lengths, and Fig. 12b showing 
the distribution of the estimation error according to this invention; 

[0036] Figures 13a and 13b illustrate the US data rate estimation results, with Fig. 13a 
showing a comparison between the actual and the estimated US date rates, and Fig. 13b 
showing the distribution of the estimation error according to this invention; 

[0037] Figures 14a and 14 b illustrate the DS data rate estimation results, with Fig. 14a 
showing a comparison between the actual and the estimated DS date rates, and Fig. 14b 
showing a distribution of the estimation error according to this invention; and 

[0038] Fig. 1 5 is a flowchart illustrating an exemplary operational flow according to this 
invention. 
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[0039] The exemplary embodiments of this invention will be described in relation to 
acquiring, forwarding, if appropriate, and analyzing diagnostic information in a 
communications environment. However, it should be appreciated, that in general, the 
systems and methods of this invention would work equally well for any type of 
communication system in any environment. 

[0040] The exemplary systems and methods of this invention will be described in relation 
to DSL modems and associated communication hardware, software and communication 
chaimels. However, to avoid unnecessarily obscuring the present invention, the following 
description omits well-known structures and devices that may be shown in block diagram 
form or otherwise summarized. 

[0041] For purposes of explanation, nxmierous details are set forth in order to provide a 
thorough understanding of the present invention, it should be appreciated however that the 
present invention may be practiced in a variety of ways beyond the specific details set forth 
herein. For example, the systems and methods of this invention can generally be applied to 
any type of communication system within any environment and for the detection of a data 
rate in any digital commimications enviroimaent. 

[0042] Furthermore, while the exemplary embodiments illustrated herein show the 
various components of the system collocated, it is to be appreciated that the various 
components of the system can be located at distant portions of a distributed network, such as 
a telecommunications network and/or the Internet, or within a dedicated secure, unsecured 
and/or encrypted system. Thus, it should be appreciated that the components of the system 
can be combined into one or more devices, such as a modern^ or collocated on a particular 
node of a distributed network, such as a telecommunications network. As wiU be appreciated 
from the following description, and for reasons of computational efficiency, the components 
of the system can be arranged at any location within a distributed network without affecting 
the operation of the system. For example, the various components can be located in a Central 
Office (CO or ATU-C) modem, a Customer Premises Modem (CPE or ATU-R), or some 
combination thereof Similarly, the functionality of the system could be distributed between 
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10043] Furthennore, it should be appreciated that the various links, includmg 
communications line 20, comiecting the elements can be wired or wireless links, or any 
combmation hereof, or any other known or later developed element(s) that is capable of 
supplying and/or communicating data to and firom the connected elements. The term module 
as used herein can refer to any known or later developed hardware, software or combination 
of hardware and software that is capable of performing the fimctionaHty associated with an 
element. 

[0044] An exemplary ranbodimcaat of the present invention focuses on estimating data 
rates for asymmetric DSL (ADSL) service. However, and in general it is to be appreciated 
that this methodology can be ^plied to estimate data dates for any digital communications 
line. 

[0045] ADSL has an upstream (US) band, within which data is transmitted ftom the CPE 
to the CO, from Tone 6 to Tone 3 1 , and a downstream (DS) band, within which data is 
transmitted from the CO to the CPE, from Tone 32 to Tone 255. The tone interval can be, for 
example, 4312.5Hz, with the ith tone corresponding to frequency (/) according to^ = 
1x43 12.5 (Hz). The TDR echo is measured by an ADSL CO modem within one frame, and is 
averaged over 10,000 frames. Each frame has 512 time samples with the sampling rate being 
2208kHz. All subscriber lines are equivalent to a single-gauge 26AWG straight loop. 

[0046] For a straight loop with a single gauge. Gauge x, let the propagation constant be 
?4, the physical length be X» then the eqmvalent equation given in Eq. (5) can be rewritten as: 

real{/^ )• X, = reali^^e)' > ^ 

yvhsK r26 and L26 is the propagation constant and loop length of the corresponding 26AWG 
equivalent loop, respectively. Let the equivalent loop length ratio between Gauge x and 
26AWG be a^o such that L^=a^'£^, then Eq. (7) becomes: 



realiy^)- a, = real{y^). (8) 
10 
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[0047] Because the propagation constant varies across frequency, a fixed ratio across 
frequency has to be computed in a least sqviare sense over a certain frequency band: 

a,{m,n)=^ , (9) 

i—m 

where^ is the frequency of the ith tone, and m and n determine the frequency band under 
consideration. The equivalent ratio changes with m and n. 

[0048] Because a subscriber line provides more attenuation on higher frequency 
components than on lower frequency components, the measured TDR echo is dominated by 
low frequency components. The spectrum of the pulse also influences the frequency band of 
the measured TDR echo. Fig. 2 shows the spectrum of the TDR pulse used in this invention. 
Fig. 3 shows the spectra of the far-end echoes from 26AWG straight loops with various 
lengths. Eighteen curves are plotted in each plot. Each cxirve corresponds to a 26AWG loop 
with a certain length. The loop length varies from Ikft to 18kft in Ikft step. It can be seen 
that the far-end echo is dominated by low frequency components, and more than 95% of the 
energy is distributed below Tone 50 (=220kHz). Since there is not much energy below Tone 
6, the TDR band is set in this exemplary embodiment from Tone 6 to 50, however, the TDR 
band can be varied as appropriate. 

[0049] As discussed above, the US band of an ADSL application is from Tone 6 to Tone 
3 1, the DS band is from Tone 32 to Tone 255, so the equivalent ratio for the US case is 
different from the DS case. Consider the difference in frequency band among the US, DS, 
and TDR case, define: 

^x(eq_us) = equivalent loop length ratio corresponding to the same US data rate, 
^xfeg_ds) = equivalent loop length ratio corresponding to the same DS data rate, and 
<^x(eqjdr) = cquivalcnt loop length ratio corresponding to the same far-end echo 
(shape and amplitude only). 
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then the equation for determining each ratio using Eq, (9) is 

axfeqjis)^ cix(m^32, n^255) , 
Clx(eqjdr)= cix(m-=6, n=^50) , 

[0050] Table 1 shows the equivalent ratios for both American loops (AWG) and 
Europeaa loops (metric). All ratios are determined using the wire primary parameters, 
characteristic impedance Zo and propagation constant y; published in ITU G. 996.1. 



Gauge X 



0.32mm 



Equivalent- 
US-Rate Ratio 



0.8 



Equivalent- 
DS-rate Ratio 



0.8 



Equivalent- 
TDR-echo 
Ratio 

Ox/ea idr) 



0.8 



Velocity 
Coefficient 
Vx (sample/^) 



6.3 



Time-shift 
Coefficient 
^ (fi/scanple) 



ANY 



0.4nun 



1.0 



1.0 



1.0 



6.9 



ANY 



O.Smin 



1.5 



1.3 



1.5 



7.5 



30 



0.63mm 



2.2 



1.6 



2.0 



6.9 



40 



0.9mm 



3.3 



2.2 



3.0 



6.5 



30 



19AWG 



3.2 



2.3 



2.8 



6.9 



20 



22AWG 



2.0 



1.6 



1.9 



7.1 



30 



24AWG 



26AWG 



1.4 



1.0 



1.27 



1.0 



1.4 



1.0 



7.1 



7.0 



Table 1: Equivalent coefficients 



40 



ANY 



[0051] Fig.4 shows the data rate comparison between 24AWG and 26AWG equivalent- 
data-rate loops. More specifically. Fig. 4a illustrates the data rate comparison between 
equivalent US data rate loops and Fig. 4b illustrates data rate comparison between equivalent 
DS data rate loops. According to Table 1, the equivalent-US-rate ratio is 1 .4; the equivalent- 
DS-rate ratio is 1.27. Fig. 4(a) shows the US data rates of 24AWG loops (length = L) versus 
that of 26AWG equivalent-US-rate loops (length = L/1 .4). Fig. 4(b) shows the DS data rates 
of 24AWG loops (length = L) versus that of 26AWG equivalent-DS-rate loops Oength = 
L/1 .27). It can be seen that the equivalent loops do have very similar data rate. The averaged 
data rate difference between equivalent loops is about 30 kbps for the US case, and 160 kbps 
for the DS case. 



[0052] Fig. 5 illustrates the far-end echo comparison between 24AWG and 26AWG 
equivalent-TDR-echo loops. According to Table 1, the equivalent-TDR-echo ratio is 1.4. 
Fig. 5a shows the comparison betweai a 26AWG 5kft loop and a 24AWG Tkft loop 
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(5* 1 ,4=7), Fig. 5b shows the comparison between a 26AWG lOkft loop and a 24AWG 14kft 
loop (1 0* 1 .4=14). It can be seen that between the equivalent loops, the amplitude and the 
shape of the far-end echoes are very similar; while the time delays are different - the 24AWG 
loop has a longer delay than the 26AWG loop. 

[0053] The equivalent ratios given in Table 1 are computed under the assumption that the 
difference in characteristic impedance Zq can be ignored. The comparisons shown in Fig. 4 
and in Fig. 5 indicate this approximation is reasonable. 

[0054] Due to the existence of gauge changes within a mixed-gauge loop, the mixed- 
gauge case is more complicated than the single gauge case. An exemplary transmission line 
system with a two-section mixed-gauge loop is schematically represented in Fig. 6. It is 
sunilar to that shown m Fig. 1 . The only difference is &e line has two sections instead of one 
section. 

[0055] Let tiie characteristic impedance, propagation constant, and loop length of the first 
section be Zqu Yi and Lu and that of the second section be Z02, n and L2. The gauge change 
point is denoted as A, and the reflection coefficient at A as p^, where: 

-^02 "j^Ol 



Zm +2^ft 



[0056] According to the transmission line theory, the wave propagation within this two- 
section loop can be represented as follows: 
Incident wave: 

1^ trip, wave at A: 

Forward wave V^^ = V^°^ -e-'^ 

Backward wave V^^ = V^""^ • e"'''^ • 



1^ trip, incident from Section 1 to Section 2: 
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1^ trip, return from the end of Section 2 to A: 

Forwardwave -e"^-^^ -(-pj-p^ 

Backward wave V^^ = V^^ ' ^'^^'^ ' Pl 

1^' trip, incident from Section 2 to Section 1 : 

1^ reflection at the source output: 

V^'' = [V^^ -e-^ +F« -e-^l-a + A) 

= fe^ • e--'^ + (V^> + Fi^-)- e-^ ]. (1 + /,J 



[0057] When each section is not too short (e.g., Li, la ^ 1000ft), higher order reflections 
can be ignored, thus the 1^ order reflection given by Eq. (10) is a reasonable approximation 
of the overall fer-end echo. Eq, 10 shows that the &r-end echo includes two dominant 
reflections, one is fiom the gauge change (the first term in Eq. (10)), the other is from the end 
of the line (the second term in Eq. (10)). Because the reflection coefficient at the gauge 
change, pA, is usually very small, /?/ is even smaller, the reflection from the end of the hue 
can be simplified to: 

ys~f:^e-^-'^-''^^'pAl^Psy (11) 



[0058] Assume the equivalent-TDR ratio of a 26AWG loop to be at(eg_uir) for Section 1 , 
and a2(eqjdr) for Section 2. Eq. (1 1) indicates that if is similar to the characteristic 
impedance of a 26AWG loop, the reflection firom the end of the mixed-gauge loop has a 
sinular shape and amplitude as the fer-end echo firom a 26AWG loop with lengtti 
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[0059] Fig. 7 shows far-end echo comparisons between a mixed-gauge loop and its 
26AWG equivalent-TDR-echo loop. More specifically, Fig. 7a shows the comparison 
between a mixed-gauge loop, which has a first section of 3kft 26AWG and a second section 
of 3kft24AWG, and a26AWG straight loop with a length of S.lkft (3+3/1.4=5.1). Fig. 7b 
shows the comparison between a mixed-gauge loop, which has a first section of 2kft 24AWG 
and a second section of 5kft 26AWG, and a 26AWG straight loop with a length of 6.4kft 
(2/1 .4+5=6.4), Fig. 7(a) suggests that ui the mixed-gauge case, although part of the far-end 
echo is the return from the gauge change, the return from the end of the loop is very sinwlar to 
the far-end echo from the equivalent straight loop. The sign of the return from the gauge 
change is the inverse of the retum from the eatid of the loop. This is consistent with tite fact 
that the gauge change is from 26AWG to 24AWG. Because the characteristic impedance of a 
24AWG loop is smaller than that of a 26AWG loop, the reflection coefficient at the gauge 
change is negative, which results in a negative retum. 

[0060] Similar to Fig. 7a, Fig. 7b also suggests that in the mixed-gauge case, although 
part of the far-end echo is the retum from the gauge change, the retum from the end of the 
loop is very similar to the far-end echo from the equivalent straight loop. However, now the 
gauge change is from 24AWG to 26AWG, the reflection coefficient at the gauge change is 
positive, so the retum from the gauge change is somewhat similar to the retum from the end 
of the loop. Therefore, it is expected that the retum fi^m tiie gauge change may have more 
influence on 24AWG+26AWG mixed-gauge case than on 26AWG+24 AWG mixed-gauge 
case. 

[0061] Assume Zs = Zoj, Zi = Z02, the voltage at the load input, Vl, is: 

therefore, the transfer fimction of the mixed-gauge line is: 

H^^^e'^^'^'^^'- -(l + pj- (12) 



15 



wo 2005/078951 



PCTAJS2005/004015 



[0062] As mentioned before, the reflection coefficient at the gauge change, pa^ is usually 
very small, thus the transfer function can be simplified to: 

H^e^^'^'^^^. (13) 

[0063] Assume the equivalent US and DS ratios to a 26AWG loops to be ajfeg_us) and 
cii(eq_d5) for Section 1, and a2(eqj4s) and a2(eg_ds) for Section 2. Eq. 13 indicates that the 2- 
section loop has the same data rate as a 26AWG straight loop with length: 

and has the same DS data rate as a 26AWG straight loop with length: 

[0064] Fig. 8 shows data rate comparison between 2-section mixed-gauge loops and their 
26AWG equivalent-data-rate loops. The first section of each loop is 26AWG and the second 
section is 24AWG. Each section varies jBx)m Ikft to 9kfl in 500fl step. It can be seen that the 
approximation given in Eq. (13) is reasonable. 

[0065] Fig. 9 illustrates the data rate comparison between 2-section loops with different 
gauge change order: one case is firom 26AWG to 24AWG; the other case is Jfrom 24AWG to 
26AWG. It can be seen that the data rates of these two cases are very similar. According to 
Eq. (12), the transfer function of the 26AWG+24AWG case be written as: 

^ 26 Am+24Am} Pa)9 

and the 24AWG+26AWG case be written as: 

^ 24Am+26Am ^ V Pa) 9 

thus the difference between these two cases is the sign of the reflection coefficient. The 
similarity in data rates between these two cases indicates the reflection at the gauge change 
has no significant influence on the data rates. Therefore, Eq. (13) is a reasonable 
approximation of the transfer function of a 2-section loop. 

[0066] Although Eq. (1 1 ) and Eq. (1 3) are derived for 2-section loops, they hold true for 
multi-section loops as weU. For a loop with n sections, assume the physical lengdi, 

16 



wo 2005/078951 



PCT/US2005/004015 



equivalent US ratio, equivalent DS ratio, and equivalent TDR ratio for the ith section (i = 1 - 
n) to be U atfeg usj, a^eq ds), and ai(eqjdrh respectively. Let the 26AWG equivalent US length, 
equivalent DS length, and equivalent TDR length be Leq_us> Leqjis» and Leqj^^ respectively, 
then: 

n 

/«! 
II 

^eq^ds 2 / ^i(eq ds) » 
n 

1=1 



[0067] As mentioned above, the ultimate goal of the equivalent loop estimation is to 
predict the data rate for a given loop. The eqmvalent lengtii estimated from the measured 
TDR echo is the equivalent TDR length. According to Table 1, the equivalent TDR ratio and 
the equivalent data rate ratio, especial in the DS case, are different, thus the data rate 
predicted using the equivalent TDR length would be inaccurate. In an exemplary 
embodiment of the present invention, the equivalent-data-rate length is derived from the 
estimated equivalent TDR length and the time delay between the measured far-end echo and 
the equivalent far-end echo. 

[0068] The time delay of a far-end echo is determined by the physical length of the loop 
and the propagation velocity of the traveling wave. Based on the transmission line theory, 
the time delay can be represented as: 

where L is the physical length of the loop, and Vp is the propagation velocity. Since a far-end 
echo is a round trip return, the numerator is double loop length. The propagation velocity of a 
transmission Ime relates to the imaginary part of the propagation constant 

V = ^ 
^ imag(y) ' 

where fi> is the radian frequency (o- 27^, and imag(j^ is the imaginary part of y. Because a 
subscriber line usually has a propagation constant whose imaginary part is a straight line 
across frequency, the propagation velocity is a constant across frequency. Thus, the 
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propagation velocity of Gauge x, denoted as can be computed fix)m the propagation 
constant as in the following 

^'^—l » (14) 

^mag\y, {f„ )]- imag\y, (/„ )]} 

where/ is the frequency of the ith tone, yx is the propagation constant of Gauge x, Af is the 
tone interval, which in the ADSL case is 43 12.5Hz, and m and « correspond to the frequency 
range of the TDR echo. In this specific case, /w = 6, w = 50. The unit for Vx is "m/s". 

[0069] The propagation velocity can be expressed in terms of number of samples per kft. 

For Gauge x, define number of samples per kft as velocity coefficient Ax. then: 

^ _ 2xl000(/() 
' 3.2808(/^/m) ' 

whereas is the sampling rate of the TDR measurement system. The unit for is "Time 
Sample/kft". The velocity coefficient of AWG loops and metric loops are Usted in Table 1 
for^ = 2208kHz. It can be seen that a 24AWG loop has ahnost the same velocity coefficient 
as a 26AWG loop, about 7 samples/kft. Consider the equivalent TDR ratio between a 
26AWG loop and a 24A WG loop is 1 .4, the delay difference, or the time shift, between a 
24AWG loop (length = L) and its 26AWG equivalent loop (length = L/L 4) is Lx7 - L/1 .4x7 
= 2xL, where Z is in kft. For the far-end echoes shown in Fig. 5 and Fig. 7, the theoretical 
time shifts between the measured far-end echo and its TDR-equivalent loop are: 

Fig. 5(a): 24AWG length = 7 kft, time shift = 14 samples; 

Fig. 5(b): 24AWG length = 10 kft, time shift = 28 samples; 

Fig. 7(a): 26AWG length = 3 kft, 24AWG length = 3 kft, time shift = 6 samples; 

Fig. 7(b): 24AWG length = 2 kft, 26AWG length = 5 kft, time shift = 4 samples. 

[0070] The time shifts given above are consistent with tiie plots given in Fig. 5 and 7. 
This indicates the time shift contains the information about the physical length of a given 
loop. 

[0071] For a straight loop with a single gauge. Gauge x, if its physical length is Z„ then 
the equivalent-TDR length Legjdr= LJax(eqjdr). Let the velocity coefficient of Gauge x be A,, 
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and that of 26AWG be /lo, then the time shift between the measured echo and the equivalent 
echo, denoted as s, is: 

where the unit of s is nimiber of samples. The physical length, Lx^ can be rewritten as: 

s 



4= 



(15) 



[0072] Eq. (15) shows that once Gauge x is known, the physical length of a straight loop, 
Zjc, can be obtained directly from the time shift s. 



[0073] For a straight loop mixed with 26AWG and Gauge x, as illustrated in Fig. 1 0, 
assume the length of the 26AWG section to be Lq and that of Gauge x to be L^^ If we know 
Xo, Ax, ax(eqjdr)y ax(eq_dsh ^® equivalcnt-TDR length Leqjdr, and the time shift s, then the 
equivalent-DS-rate length, denoted as Leqjis^ can be obtained by solving the following 
equations: 

Leq_tdr = A + ^x/ ^x{eq^tdr) 
' ^eq^ds = ^0 + ^x/ ^x{eq^ds) * 
S " VV'^' ""^o/^xCeg.lcfr)) 



[0074] The solution is: 



K'Kl^xieqJdr) \j^x{eq^ds) ^x{eq_idr) _ 



(16) 



[0075] Define the time-shift coefBcient for Gauge x be ^, i.e. 



1 



1 



K^^/^x{eq^tdr) \j^x{eq^ds) ^x(eg^idr) ^ 



then Eq. (16) can be written as: 



Leg_ds = ^eg^tdr "^^'^x' 



(17) 



[0076] Eq. (17) shows that the equivalent-data-rate length, Leqjis, can be determined 
based on the equivalent-TDR-echo length, Leqjdn the time shift s, and the time shift 
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coefficient ^. Table 1 lists the time shift coefficients of both AWG wires and metric wires. 
Because for 0.4mm and 26AWG wires, the time shift is 0 for any loop length, in other words, 
the equivalent-data length, the equivalent-TDR-echo length and the physical length are the 
same for these wires, the time-shift coefficient for these two gauges can be any value. Eq. 
(17) indicates that if the time shift coefficient, z^c, is a constant across all gauges, the 
relationship given by Eq. (17) would be independent of Gauge jc. However, Table 1 shows 
the time-shift coefficients are not identical across gauge. Since the most popiilar gauges used 
in the field are 24AWG and 26AWG in North American, and 0.4mm, 0.5mm, and 0.63mm in 
Europe, we average the time shift coefficient across 24AWG, 0.5mm and 0.63mm, the 
rounded average is i^ftean ^ 40. 

[0077] Because Ihe equivalent-US-rate ratio is very similar to the equivalent-TDR ratio, 
no correction is made for US rate prediction. 

[0078] As discussed, the goal of the equivalent estimation method is to predict the data 
rate for a given subscriber line based on a TDR measurement. The input of the method is a 
measured TDR echo, the output is the predicted DS and US data rates. The intermediate 
steps include the equivalent TDR length estimation and the length correction for data rate 
prediction. In order to predict the data rate correctly, this exemplary method needs to know 
the data rate versiis loop length curve of 26 AWG straight loops. The detailed procedure is 
discussed below in relation to Figs. 1 1 and 12. 

[0079] Fig. 1 1 illustrates an exemplary channel capacity estimator 100 according to this 
invention. The estimator 100 comprises a TDR measurement module 1 10, a theoretical TDR 
echo determination modide, an equivalent TDR length estimation modulel30, an optimizer 
module 140, a memory 150, a controller 160, a correction module 170 and an US and DS 
data rate prediction module 180, all interconnected by link 5. The estimator 100 is connected 
to a line 20, which is in turn connected to a modem 200, such as a DSL modem. 

[0080] In operation, the TDR measxirement module 110 reads a previously measured or 
detemiines the measured TDR echo, (echojneasured(i)\ where / is the time sample index. 
In the case of ADSL, one fi:ame has 512 samples, thus z = 0 - 511. 
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[0081] Next, the theoretical TDR echo deteimination module 120 determines the 
theoretical TDR echoes for 26AWG straight loops with various loop lengths. 

[0082] Let the loop length of the nth loop be X„, w = 7 - where N is total number of 
loops, the corresponding theoretical TDR echo be echojnodeln(i), then: 

" Zo/tanhQ^:,) 
_Z,+Zo/tanh(/Zj 

where Zq and /are the characteristic impedance and propagation constant of a 26AWG loop, 
Zs and ^ are the output impedance and voltage of the source. 

[0083] Then, the equivalent TDR length estimation modide 130 estimates the equivalent 
TDR length. This is achieved by: 

(1) Finding the best time shift between the measured far-end echo and the theoretical 
far-end echo of the nth loop by solving the following optimization problem: 

min ^ {echo _ measured{j + m)- echo _ model„ {j)f , 

"* y=o 

where w is a variable representing the time shift, which varies in a certaia region 
with a one-sample step. The minimum error across m is denoted as E(nX and the 
corresponding best shift as S(n). 

(2) Let n = w+i. If w <N, go to (1); otherwise go to (3). 

(3) Let n* = min£(w), then the equivalent TDR length. Leg tdr, is: 

n — 

and the corresponding time shift, denoted as j*, is: 

[0084] The correction module 170 then corrects the equivalent-TDR length for data rate 
predictioxL 

[0085] According to the analysis given above, the equivalent-US-rate length and the 
equivalent-DS-rate length relate to Leqjdr and j?*: 



21 



wo 2005/078951 



PCT/US2005/004015 



*eqJLdr > 



and 



'eqjdr 



mean 



'eqjtdr 



+ 5 * AQ(Ji I sample) . 



[0086] The US and DS data rate prediction module 1 80 then predicts the US and DS data 
rates, which can be one or more of output and or displayed on a display device (not shown), 
by letting the US rate-length function for 26AWG loops be Rateus(L), and the DS rate-length 
function be RateDs(L), then the US rate, denoted as US^Rate^ and the DS rate, denoted as 
DS_Rate^ are: 



[0087] The present invention has been tested on several different CO modems that have 
TDR functionality. The results from one of the modems are given below. All of the loops 
tested are listed in Table 2. The total number of loops is 1291 . All loops are straight loops 
with a single section or up to four sections. Both American wires and European wires are 
tested. Fig. 12 shows the estimation results on equivalent-TDR length, with Fig. 12a showiag 
the estimated length versus theoretical length, which is calculated using Table 1, and Fig. 12b 
showiag the distribution of the estimation error. 

[0088] Fig. 13 shows the estimation results on US data rate, with Fig. 13a showing the 
estimated US rate versus the measured data rate, i.e. the actual data rate when connecting a 
CO-CPE modem pair using the given loop, and Fig. 13b showing the distribution of the error 
on US rate estimation. 

[0089] Fig. 14 shows the estimation results on DS data rate, with Fig. 14a showing the 
estimated DS rate versus the measured data rate, i.e. the actual data rate \^en connecting a 
CO-CPE modem pair using the given loop, and Fig. 14b showing the distribution of the error 
on DS rate estimation. 



US_Rate = Rate^ (x,^^ ) , and 
DS^Rate = Rate^ (x^^^ ) . 
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No. of 
sections 


Configuration 


Length 


No. 
of 

loops 


1 -section 


24A WG (L) 


L = Ikit - 22kft / Ikft 


22 


26A WG (L) 


L= Ikft- 18kft/ 100ft 


171 


0.4 mm (L) 


L = 200m - 3500m / 50m 


67 


0.5 mm (L) 


L = 200m - 5000m / 50m 


97 


2-section 


26AWG (LI) + 24AWG(L2) 


LI, L2 = Ikft - 9kft / 500ft 


289 


0.4mm(Ll) + 0.5mmOL2) 


LI, L2 = 200m - 3000m / 200m 


225 


ETSI #3 (L=length of 0.4mm 
section) 


L = 50m - 3000m / 50 m 


60 


3-section 


26awg(Ll)+24awg(L2)+26awg(L3 
) 


LI, L2, L3 = Ikft - 9kft / 2kft 
Ll+L2+L3<=18kft 


90 


ETSI #6 (L^length of 0.4mm 
section) 


L = 50m - 3000m/ 50 m 


60 


ETSI #7 (L=length of 0.4mm 
section) 


L = 50m - 3000m / 50 m 


60 


4-section 


26awg+24awg+26awg+24awg 


The length of each section is 
ranaomly Chosen. Ine total 
physical length is less than 18kft. 


30 


ETSI #4 (L=length of 0.4mm 
section) 


L = 50m - 3000m /50 m 


60 


ETSI #5 (l>length of 0.4mm 
section) 


L = 50m -3000m /50 m 


60 



Table 2: A list of exemplary tested loops 



[0090] Exemplary Figs. 12, 13 and 14 indicate that the equivalent-TDR length estimation 
error is less than 500ft for 96% of the loops; that the US data rate estimation error is less than 
100kbps for 99% of the loops, and that the DS data rate estimation error is less than 500kbps 
for 97% of the loops. 

[0091] Fig. 15 illustrates an exemplary method of predicting data rates according to this 
invention. In particular, control begins in step SI 00 and continues to step S200. In step 
S200, a measured TDR echo, is obtained. Next, in step S300, tiie theoretical TDR echo for 
26AWG straight loops with various loop lengtiis are determined and stored. It is to be 
appreciated howevCT that the gauge of the equivalent loop does not necessary have to be 
26AWG, to the contrary, the gauge could be any gaiige. The only requirement is the data rate 
of the selected gauge should be known. Furthermore, the described methodology is not only 
applicable to single gauge straight loops, but also multi-section loops with different gauges as 
well as loops wifli bridged taps. 
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[0092] Then, in step S400, the equivalent TDR length is estimated. As discussed above, 
this is detennined by: 

(1) Finding the best time shift between the measured far-end echo and the 
theoretical far-end echo of the nth loop by solving the following optimization 
problem: 

WN-A 

min ^ \echo _ meastired(j + m)- echo _ model „ {j)f , 

y=o 

(2) where mis a variable representing the time shift, which varies in a certain 
region with a one-sample step. The minimxmi error across m is denoted as 
E(n), and tiie corresponding best shift as S(nX 

(3) Let n = n+l. If n <N, go to (1); otherwise go to (3). 

(4) Let «* = min£:(«), then the equivalent TDR length. Leg tdr, is: 

n — 

(5) and the corresponding time shift, denoted as s^^ is: 
Control then continues to step S500. 

[0093] In step S500, the equivalent-TDR length for data rate prediction is corrected. 
According to the analysis given above, the equivalent-US-rate length and the equivalent-DS- 
rate length relate to Leqjdr and 5'*: 

^eg^us " ^eqjidr » 
J^cg^ds = L^^tdr + * '-^mean = L^^tdr + * ^40(Jt / SCmpIc) . 

[0094] The US and DS data rates are then predicted in steps S600 and S700, respectively, 
in accordance with the following, which can be one or more of ou^ut and or displayed on a 
display device (not shown), by letting the US rate-length fimction for 26AWG loops be 
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Rateus(L), and the DS rate-length function be RateDs(L), then the US rate, denoted as 
USJRate^ and the DS rate, denoted as DSJtate, are: 



One or more of the US and DS data rate(s) are then output in step S800 and control continues 
to step S900 where the control sequence ends, 

[0095] The above-described system can be implemented on wired and/or wireless 
telecommunications devices, such a modem, a multicarrier modem, a DSL modem, an ADSL 
modem, an XDSL modem, a VDSL modem, a linecard, test equipment, a multicarrier 
transceiver, a wired and/or wireless wideAocal area network system, a satellite 
commimication system, a modem equipped with diagnostic capabilities, or the like, or on a 
separate programmed general purpose computer having a commxmications device. 

[0096] Additionally, the systems, methods and protocols of this invention can be 
implemented on a special purpose computer, a programmed microprocessor or 
microcontroller and peripheral integrated circuit element(s), an ASIC or other integrated 
circuit, a digital signal processor, a hard-wired electronic or logic circuit such as discrete 
element circuit, a programmable logic device such as PLD, PLA, FPGA, PAL, modem, 
transmitter/receiver, or the like. In general, any device capable of implementing a state 
machine that is in turn capable of implementing the methodology illustrated herein can be 
used to implement the various communication methods, protocols and techniques according 
to this invention. 

[0097] Furthermore, the disclosed methods may be readily implemented in software 
using object or object-oriented software development environments that provide portable 
source code that can be used on a variety of computer or workstation platforms. 
Alternatively, the disclosed system may be implemented partially or fiilly in hardware using 
standard logic circuits or VLSI design. Whether software or hardware is used to implement 
the systems in accordance with this invention is dependent on the speed and/or efQciency 
requirements of the system, the particular ftmction, and the particular software or hardware 




DS^Rate = Raters (^eg^a. ) • 
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systems or microprocessor or microcomputer systems being utilized. The communication 
systems, methods and protocols illustrated herein however can be readily implemented in 
hardware and/or software using any known or later developed systems or structures, devices 
and/or software by those of ordinary skill in the applicable art from the functional description 
provided herein and with a general basic knowledge of the computer and telecommunications 
arts. 

[00981 Moreover, the disclosed methods may be readily implemented in software, that 
can be stored on a storage medium, executed on programmed general-purpose computer, a 
special purpose computer, a microprocessor, or the like. In these instances, the systems and 
methods of this invention can be implemented as program embedded on personal computer 
such as JAVA® or CGI script, as a resource residing on a server or computer workstation, as 
a routine embedded in a dedicated communication system or system component, or the like. 
The system can also be implemented by physically incorporating the system and/or method 
into a software and/or hardware system, such as the hardware and software systems of a 
communications transceiver and operations support system. 

[0099] It is therefore apparent that there has been provided, in accordance with the 
present invention, systems and methods for estimating channel data rate. While this 
invention has been described in conjimction with a number of embodiments, it is evident that 
many altematives, modifications and variations would be or are apparent to those of ordinary 
skill in the applicable arts. Accordingly, it is intended to embrace all such altematives, 
modifications, equivalents and variations that are within the spirit and scope of this invention. 
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